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Abstract 

The temperature induced phase transition is investigated in the 
one-component scalar field cj) 4 model on the lattice. Using the GPU 
cluster a huge amount of Monte Carlo simulation data is collected for 
a wide interval of coupling values. This gives a possibility to deter- 
mine the low bound on the coupling constant Ao when the transition 
happens and investigate its type. We found that for the values of A 
close to this bound a weak-first-order phase transition takes place. It 
converts into a second order one with the increase of A. A compari- 
son with the results obtained in analytic and numeric calculations by 
other authors is given. 
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1 Introduction 

Scalar field models with a spontaneous symmetry breaking (SSB) are con- 
sidered in various fields of physics, like quantum field theory, collective phe- 
nomena, quantum dots, high-temperature superconductivity, etc. They often 
serve as toy models to develop both the analytic calculation schemes and the 
numeric simulation techniques to describe a wide class of temperature in- 
duced phase transitions. Mult i- component scalar field with the orthogonal 
symmetry O(N) is a popular choice for these investigations. The simplest 
model contains one-component scalar field and can be called the 0(l)-model. 

There is a long history of studying phase transitions in the O(N) models 
(see Refs. p]-[3] and references therein). In analytic calculations based on the 



1 



perturbation theory (PT) methods, various resummation schemes are used. 
However, the application of different resummation techniques leads to con- 
tradictory results about the type of the phase transition. In Ref. [5] a second 
order phase transition for the 0(iV)-model was determined independently of 
the coupling value by applying some kind of resummations. The same result 
was also derived by using the renormalization group approach [B]. On the 
contrary, for the case of 0(l)-model the phase transition of first order was 
observed in the daisy, super daisy and some type beyond resummations for 
the extremely weak coupling constant [I]. It was shown that various kind 
resummations can lack their expansion parameters near the phase transition 
temperature T ~ T c . The first order phase transition was also observed 
within the 2PI formalism in the double-bubble approximation Ref. [7]. 

The discrepancies in the analytic results can be resolved by applying the 
Monte Carlo (MC) simulations on the lattice. As numerous MC simulations 
showed, the second order phase transition takes place. So, nowadays a gen- 
eral belief is that the phase transition is of second order and PT fails in 
this problem. However, all the performed already MC calculations used the 
coupling values A > 0.01. A weak- first-order phase transition determined in 
Ref. [1] for an extremely weak couplings, has never been determined on the 
lattice. This notion - 'extremely weak coupling' - assumes some physical mo- 
tivation, namely the so-called Linde- Weinberg bound on the scalar field mass 
[S]) [IH]- For many years ago these authors have observed that in the mod- 
els with negative mass squared m? < the SSB does not happen for small 
values of the coupling constant A < Ao even at zero temperature. Although 
the actual value Ao depends on the mass parameter entering Lagrangian, 
it is natural to consider small coupling values to reach the Linde- Weinberg 
bound. Physically, this bound reflects an important property of the SSB - 
the existence of the range of parameters allowing the total effective poten- 
tial to be dominated by the positive radiation quantum effects instead of 
the negative classical part. In this sense, the notion 'Linde- Weinberg bound' 
is reasonable independently of a particular model considered. In fact, the 
'Linde- Weinberg bound' A was never discussed within MC simulations, so 
the 'small' values A ~ 0.01 commonly used in the simulations can occur much 
greater than Ao- The goal of the present paper is to fill the gap and to an- 
swer the question about the type of the phase transition in the 0(l)-model 
in MC simulations for values of the coupling constant much below the values 
investigated already. 

A well known approach to determine the type of a phase transition is to 
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investigate the behavior of the order parameter. For the first order transition 
overheated and supercooled meta-stable states arise at the critical tempera- 
ture dependently of the initial configuration. With the initial ordered phase 
(cold start) the overheated configurations dominate. In contrast, with the ini- 
tial disordered phase (hot start) the supercooled configurations are mainly 
observed. Since the order parameter distinguishes these alternative meta- 
stable phases, collecting the statistics for different starts together one can 
see a hysteresis plot for the first order transition. Of course, no hysteresis 
can be found for the order parameter in case of the second order transition. 
These technique was successfully applied, for example, to determine the order 
of the phase transition in the lattice QCD [TTj . 

In case of the 0(l)-model, the evident order parameter is the field con- 
densate taking non-zero values in the phase with the broken symmetry. In 
MC simulations the condensate can be easily measured as the field average 
over the lattice. In the present paper, we collect the statistics for a wide 
interval of coupling values searching for possible hysteresis behavior. As a 
result, at A > 10~ 3 no hysteresis is observed confirming the second order 
phase transition. However, at A < 10~ 3 we find the hysteresis behavior, 
and the hysteresis becomes more pronounced with A decreasing. Thus, we 
conclude that the type of the phase transition changes at extremely weak 
couplings. Moreover, at A ~ 10~ 5 the ordered phase does not occur with the 
hot start even at zero temperature indicating, probably, the Linde- Weinberg 
bound A . 

The paper is organized as follows. In Sect. 2 we introduce a parametriza- 
tion of the 0(l)-model on a lattice allowing to produce stable results near 
the critical temperature in a wide interval of couplings. Sect. 3 contains the 
results of MC simulations. Sect. 4 is devoted to conclusions. 



In continuous space the thermodynamical properties of the model are de- 
scribed by the generating functional 



2 The model 





where if is real scalar field, and the action is 
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The standard realization of generating functional in Monte Carlo simulations 
on a lattice assumes a space-time discretization and the probing random val- 
ues of fields in order to construct a Boltzmann ensemble of field configura- 
tions. Then any macroscopic observable can be measured by averaging the 
corresponding microscopic quantity over this ensemble. 

The 0(l)-model on a lattice contains three energy scales: the mass m, the 
temperature T, and the inverse lattice spacing. The mass m can be chosen 
as a unit, then the temperature (together with A) determines the physically 
important values of ip, and the lattice spacing must support successful sim- 
ulation of these field values. Since the field p is distributed in the infinite 
interval (—00,00), the characteristic scale can also move in an extremely 
wide interval when A changes. In principle, MC algorithms are able to find 
an unknown scale in the infinite interval. However, such a strategy requires 
some non-trivial parameter re-setting for different values of A. Instead, we 
prefer to rewrite the model in terms of dimensionless variables taking values 
from finite intervals, automating also the correct choice of the lattice spacing. 

First we introduce one-to-one transformation tp(U) to new field variable 
U(x) defined in the finite interval (0, 1). Let U = 0.5 corresponds to ip — 
and <p(U) = — ip(l — U). The generating functional in terms of U reads 



For Monte Carlo simulations we introduce a hypercubic lattice with hyper- 
torous geometry. We use an anisotropic cubic lattice with a spatial and a 
temporal lattice spacing a s and at = a s /( with £ > 1, respectively. The 
scalar field is defined in the lattice sites. Transforming the Jacobian as 
det A = exp (Tr log A), the generating functional becomes 
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The lattice forward derivative is defined as usually by the finite difference 
operation 



where a M is the lattice spacing in the \i direction, fi is the unit vector in the 
direction indicated by /i. 

In the case of pure condensate field the action is determined by the po- 
tential 



V[U] 



(7) 



This potential is topologically equivalent to the potential V(cp) = —m 2 {p 2 /2 
+X(p 4 /4:. It has one local maximum at U — 0.5 and two symmetric global 
minima at Uq and 1 — Uq. The spread between the values of the potential at 
local maximum and global minima is 

The quantities Uq and Ay play a crucial role in Monte Carlo simulations. 

Considering the phase transition, one must guarantee that the Monte 
Carlo algorithm meets the field values compatible with both the phases to 
choose. If Uq — » 0.5, then the broken phase can be missed since the corre- 
sponding field values are extremely rare events. On the other hand, Uq — > 
{Uq — »■ 1) washes the unbroken phase out. So, to study the phase transition 
in the model, we choose the following conditions: 

Uq = 0.25, Ay = 1. (9) 

Thus, the half of generated field values will be between the global minima 
of the 'effective' potential, and no phase will be accidentally missed. The 
probability to prefer condensate or non-condensate values will be of order 
~ 0.5 ensuring the fast convergence of Monte Carlo algorithm. As it will 
be shown, these conditions successfully works for 10 -5 < A < 0.5. For the 
larger (smaller) A the condensate becomes too weak (strong) and the numeric 
values in ([9]) must be reconsidered. 
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To satisfy two conditions we use a convenient two-parameter function 



<p[U\ 



m£arctanh 7/(2*7 - 1) + (1 - 77) (2*7 - 1) 



(10) 



with £ > and 77 > 0. The values of £ and 77 have to be found as the solution 
of equations dV / dU\u=u = and (jSJ). These equations can be written as 



2F[u ] - g[u ] 

(F[u ]-g[u Q }f 
x 



AC 
m 4 af 



m 2 z 



<p2[u ] = jr[u ] - g[u ] 



where z is a dimensionless parameter of the model, 

1-7/(2*7-1) - (1 -77) (2*7- l) 3 



K[U\ 
F[U) 



l + r/(2f/-l) 
( K"[U\K,[U] _ 
V (/C'p7])» 
-/C'[*7] 



f (l-77)(2*7 
1] log/C[*7], 



4 log 



4t7/C[*7] 



At 
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(13) 
(14) 
(15) 



where the primes denote derivatives. Eq. ( flTT) gives r)(z), then £ can be found 
from ffl2|) . There is no physical solution for 2; < z m i n . This forbidden interval 
corresponds to low temperatures which cannot be reached within the chosen 
parametrization. Finally the lattice action is 



S[U(x)} 
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*7(x + a M /t) — *7(x) 



--£[*7(x)] - Flog 2 /C[*7(x)] + F 2 ,zlog 4 /C[*7(x)] + V 
\-GPo] 



21og z /C[*7 



-Ay - log(2m£77). 



(16) 
(17) 



The constant part of the action Vq is completely unimportant for calculations 
and can be omitted, since Monte Carlo algorithm is based on the difference 
between the actions of modified and initial field configurations. 

By varying ( it is possible to change a t , while keeping a s fixed. Conse- 
quently the temperature T ~ ( can be changed continuously at fixed a s . The 



6 



2 1) 



s 



Figure 1: \<p\ in the units of classical condensate m/y/X for A = 5 ■ 10 4 on 
lattice 16 4 

field condensate, (p, is measured as the average of ip(x) over the lattice. In 
Fig. [T] we plot \tp\ in the units of classical condensate mj v^A for A = 5 ■ 10~ 4 
and 16 4 lattice. Lower values of ( and z corresponds to lower temperatures. 
One can see evident phase transition with the field condensate growing with 
the decreasing temperature. Since the clear positive or negative values of (p 
appear in lattice configurations in the broken phase, we conclude about the 
absence of domains and will use \(p\ in plots. 

To determine the type of the phase transition we consider in details three 
slices of two-dimensional function |^|(z,£) at fixed z = 0.35, z = 0.5 (for 
A = 5 • 10~ 5 ), and £ = 1. For example, the slices are shown in Fig. [I] as bold 
red lines (A = 5 • 10~ 4 ). We compute the field condensate with the hot and 
cold starts for different (. A hysteresis behavior means a first order phase 
transition. 

3 The Monte Carlo simulation results 

To estimate the order of the phase transition a large amount of simulation 
data must be prepared. The simulation requires a fairly powerful computing 
resources, especially for large lattices. To speed up essentially the generation 
of data we apply a GPU cluster as a computational platform. It consists 
of ATI Radeon GPUs: HD6970, HD5870, HD5850 and HD4870 with the 
peak performance up to 11 Tflops. The low- level AMD Intermediate Lan- 
guage (AMD IL) is used in order to obtain the maximal performance of used 
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Figure 2: Temperature dependence of the field condensate \(p\ for the lattices 
16 4 and 32 4 at z = 0.35 and z = 0.5 for C = [1.5; 2.5]. 

hardware. A trivial parallelization scheme is implemented for cluster com- 
putation. Some technical details of MC simulations on the ATI GPUs and a 
review of the AMD Stream SDK are given in Ref. p2] and references therein. 

In the MC simulations, we use lattices of different sizes up to 64 4 . Most 
statistics are obtained for lattices 32 4 and 16 4 , the qualitative behavior are 
checked on larger lattices. RANLUX pseudo-random number generator is used 
in the MC kernel, all the key results are checked with RAN MAR generator 
[15] . Lattice data are stored with a single precision. MC updating are also 
performed with the single precision whereas all averaging measurements are 
performed with the double precision to avoid the accumulation of errors. 

The system passed 5000 MC iterations for every run to be thermalized, 
then we used 1024 MC configurations (separated by 10 updates) for measur- 
ing. 



The temperature dependence of \<p\ for the lattices 32 4 and 16 4 at z = 0.35 
and z = 0.5 for £ = [1.5; 2.5] is present in Fig. |2j The whole data set for 
every plot is divided into 15 bins. Different initial conditions are shown as: 

• hot start - red circles (lattice 32 4 ) and brown dashes (lattice 16 4 ); 



S 



', a - lattice: 32" (cold start) 

! o - lattice: 32' (hot start) ~ 
01 ; x - lattice: 16' (cold start) 
- - lattice: 16' (hot start) 

00 L, 1 1 1 1 1 1 1 1 1 1 1 1 1 iM-i l^EJ H- i 

0.35 0.36 0.37 0.38 Z 

Figure 3: \<p\ in the units of classical condensate m/\/\ for A = 5 • 10 -4 and 
( = 1 on lattices 16 4 and 32 4 

• cold start - blue triangles (lattice 32 4 ) and green crosses (lattice 16 4 ). 

The mean values and 95% confidence intervals are presented by corresponding 
pointers for each bin. Every bin contains 150 points. 

As it is seen from Fig. [2J for A = 0.01 the temperature dependence of 
field condensate is not sensitive to the initial configuration. Both cold and 
hot starts lead to the same behavior of field condensate for various (. This 
meas a second order phase transition. 

For smaller values of A the overheated configurations occur in the broken 
phase for the hot start. That is, different starts demonstrate a hysteresis 
behavior. This corresponds to a phase transition of the first order. 

With further decreasing of A, for A < Ao ~ 10~ 5 the behavior of cold and 
hot starts is completely separated and independent of the temperature. This 
is plotted in Fig. [2j Such type property probably means that the SSB does 
not happen even at zero temperature and corresponding value of A can be 
identified with the Linde- Weinberg low bound. 

As it is seen in Fig. [21 two different slices z = 0.35 and z = 0.5 demon- 
strate the same behavior. The calculations for the slice ( = 1 reproduce the 
described above results again (see Fig. [3]). Thus, the hysteresis occurs at 
small A independently of the model parameters z and (. 

4 Conclusion 

As it was discovered in the MC simulations, the temperature phase transition 
in 0(1) 4 model is strongly dependent on the coupling value A. There is a low 
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bound A = determining the range where SSB is not realized. Close to 
this value in the interval 10~ 5 < A < 10~ 3 the phase transition is first order. 
For larger values of A the second order phase transition happens. These 
types of the behavior have been determined on the lattices of different sizes 
independently of the internal model parameters. Our calculation procedure 
was developed to accelerate the MC procedure in the domain of parameters 
close to the transition for a wide range of coupling. For usually considered 
values of A ~ 0.01 — 0.1 it gives the results coinciding with the ones existing 
in the literature and signalling the second order phase transition. To our 
knowledge, systematic investigations for smaller values of coupling have not 
been carried out yet. 

Our observations, in particular, may serve as a guide for the applicability 
of different kind resummations in perturbation theory. In fact, we see that 
the daisy and super daisy resummations give qualitatively correct results for 
small values of A. For larger values they become non-adequate to the second 
order nature of the phase transition. In this case other more complicated 
resummation schemes should be used. 

The change of the phase transition type dependently of the coupling value 
is not a new phenomenon. For instance, in the standard model of elementary 
particles it is well known that the electroweak phase transition is of first 
order for small A and it converts into a cross-over or even second order one 
for sufficiently large values of A. We have observed that this happens even 
in the simple model with one coupling. 

In the present investigation, we concentrated mainly on the qualitative 
aspects of converting the phase transition type due to the change of the 
coupling values. So, we skip an ubiquitous procedure relating the lattice and 
physical variables, as unessential. 
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